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Abstract 

We study the joint exit probabilities of particles in the totally asymmetric simple exclusion 
process (TASEP) from space-time sets of given form. We extend previous results on the space-time 
correlation functions of the TASEP, which correspond to exits from the sets bounded by straight 
vertical or horizontal lines. In particular, our approach allows us to remove ordering of time 
moments used in previous studies so that only a natural space-like ordering of particle coordinates 
remains. We consider sequences of general staircase-like boundaries going from the northeast to 
southwest in the space-time plane. The exit probabilities from the given sets are derived in the 
form of Fredholm determinant defined on the boundaries of the sets. In the scaling limit, the 
staircase-like boundaries are treated as approximations of continuous differentiable curves. The 
exit probabilities with respect to points of these curves belonging to arbitrary space-like path are 
shown to converge to the universal Airy2 process. 
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I. INTRODUCTION 



Consider the system of particles on the ID integer lattice. At any time moment a con- 
figuration of particles is specified by a set of N strictly increasing integers, (x\ > x 2 > ■ ■ ■ ), 
denoting particle coordinates. They evolve in a discrete time t G 7L according to the TASEP 
PQ dynamical rules: 

I. A particle takes a step forward, (scj — > 2^ + 1), with probability p and stays at the 
same site, (xj — > Xi), with probability q = 1 —p provided that the target site is empty, 

(Xi + 1 ^ JCi_i). 

II. If the next site is occupied, (xi + 1 = the particle stays with probability 1. 

III. The backward sequential update is used [2]: at each time step the positions Xi of 
all particles are updated one by one, in the order of increasing of particle index: 
z = l,2,3,... 

These dynamical rules define transition probabilities for a Markov chain constructed on 
the set of particle configurations. Given initial conditions, one can inquire for probabilities of 
different events in course of the Markov evolution. In present paper, we are interested in the 
correlation functions which are the probabilities for events associated with a few specified 
particles and given space-time positions. 

A. Spacial correlation functions of the TASEP. 

The first exact result on correlation functions in TASEP goes back to prominent Johans- 
son's work j3], where he considered the evolution of TASEP with parallel update and step 
initial conditions, 

x° = Z< , (LI) 

and obtained the distribution, Pt{xN > M — N), of the distance M traveled by N-th particle 
up to time t. This result was later generalized to the backward sequential update [1] and the 
flat initial conditions [5]. The connection of the TASEP with the theory of determinantal 
point processes revealed in [6], [7] allowed also calculation of the multi-particle correlation 
functions, i.e. distribution P t (x ni > ai, . . . x nm > a m ), of positions of m selected particles at 
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fixed time t, where 1 < n\ < ■ ■ ■ < n m are m integers numbering the selected particles. The 
multi-particle correlation functions were extensively studied for different initial conditions 
in a series of papers [7HH] • The result can generally be represented in a form of the Fredholm 
determinant of the operator with some integral kernel. An asymptotic analysis of the kernel 
is of special interest as it allows one to study the scaling limit of the correlation functions, 
which is believed to yield universal scaling functions of the Kardar-Parisi-Zhang (KPZ) 
universality class [TO] . 

There is a law of large numbers, which implies that the stochastic evolution converges 
to a deterministic limit [Uj [T2] . Specifically, in the TASEP, if we measure coordinate x n of 
n-th particle at time t, the deterministic relation between rescaled variables 

v = n/L,u) = tj L, ,7 = {x n + n)J L (1.2) 

holds with probability one as L — > oo. An explicit form of this relation can be found from 
the hydrodynamic conservation law 

d t p + d xJ = 0. (1.3) 

for the density of particles p. Here j = j(p) is the stationary current of particles, which is a 
model-dependent function of the density. In the case of backward update the current is 

= 

1 -pp 



Then, the solution of (1.3) with initial conditions (1.1), yields relation 



puj - ^fqv - y/j = 0, (1.5) 



which holds in the range — p/q < (7 — v)ju) < p. For the the formula (1.4) and its relation 



to (1.5) we address the reader to references [3111]. 

An exact calculation of the correlation functions allows one to study fluctuations of the 
random variables near their value on the deterministic scale. Given v and u, let 7(0;, u) be 
the rescaled particle coordinate. The deviation 5x n = x n — L(j(co, v) — u) of the particle 
coordinate x n develops on the KPZ characteristic scale fluctuations 

5x n ~ L a ,a = 1/3. (1.6) 

The distribution of the rescaled variable 

s = k' 1 lim 5x n L~ a , (1.7) 

L— >oo 



is a universal scaling function of the KPZ class, dependent only on the form of the initial 
macroscopic density profile. Note that the model dependence is incorporated into a single 
non-universal constant k x . The examples of distributions obtained from the asymptotic 
analysis of the one-point correlation function are the Tracy- Widom functions F\ and F2 for 
flat and step initial conditions respectively. These functions are well known for appearing in 
the theory of random matrices as the distributions of the largest eigenvalue in the orthog- 
onal and unitary Gaussian ensembles [131 EH] • Their presence turns out to be a universal 
feature of the KPZ class. Furthermore, the study of multipoint distributions shows that 
the fluctuations of coordinates of different particles, say x ni and x m , remain non-trivially 
correlated random variables on the scale 

|ni — n 2 \ ~ LP, (3 = 2/3. (1.8) 

This is the second power law characterizing the KPZ class. The critical exponents a = 1/3 
and (3 = 2/3 are called fluctuation and correlation exponents respectively. After correspond- 
ing rescaling of particle numbers, one arrives at the one-parametric family of correlated 
random variables: 

s{u) = hm — , (1.9) 



where K n is another non-universal constant. For the cases of flat and step initial conditions, 
the joint distributions of these variables define universal Airyi [6] and Airy2 [15] ensembles, 
whose one-point distributions are F\ and F 2 . 



B. Space-Time correlations and mapping to the last passage percolation. 

So far we have been discussing only the spacial correlations between positions of different 
particles at a fixed time moment. However, generally, one can consider joint probability 
distributions of events associated with different particles, positions and time moments, which 
happen in course of the TASEP evolution. We will refer to these distributions as the space- 
time correlation functions. An example of such a function, the distribution of positions 
of a tagged particle at different moments of time, has been calculated in [16]. A more 
general correlation function, the distribution P(x ni (ti) > ai, . . . , x nm {t m ) > a m ) of positions 
x ni , . . . , x nm of selected particles with numbers 

rii < ■ ■ ■ < n m (1.10) 
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at time moments ti, . . . , t m , was studied in [TTJIIH]. The method was used that restricted the 
analysis to the sets of space-time points, such that the time coordinates decreased weakly 
with the particle number and vice versa: 



U > U +1 , if rii < n i+1 , (1. 11) 

U > t i+1 , if Hi = n i+1 . (1.12) 

This arrangement of time moments was named space-like by the authors of [TTIIT8] . Another 
example of the space-time correlation function, the current correlation function, was recently 
obtained in [19J. This was the probability distribution P(t ni < a 1; . . . ,t nm < a m ) of time 
moments t ni , . . . , t nm at which m selected particles with numbers 

ni < • • • < n m (1.13) 

jump from the respective sites x ni , . . . , i„ m selected from the set 

{Xi = x - i + N : i = 1, . . . , n m }, (1.14) 

given x G N > n m and the initial configuration if = 1 - i,i G N. Due to non-crossing of 
space-time particle trajectories, the range of time moments accessible for the dynamics is 

t m <---<t nm . (1.15) 



The time orderings ( I.ll|1.12 ) and ( I.13[1.15 ) are opposite to each other. These orderings, 



however, have different origins. In [16HT8], numbers of particles n\, ri2, ■ ■ ■ and time moments 
t ni ,t n2 , . . . are fixed, and particle coordinates x ni , x n2 , . . . are random variables. In the case 
of current correlations [19] , time moments t ni are random, while particle coordinates x ni and 
numbers n» are related fixed parameters. Therefore, unlike ( I.ll|1.12 ) in [TTj [18]. (1.15) from 



[19] is not an external constraint, but is the consequence of dynamics: it shows domains 
which can be reached in the random process with nonzero probability. 

Which variable is chosen to be random is, however, not important in the scaling limit, 
when the three variables, time and space coordinate and the number of a particle, acquire 
equivalent significance due to separation of fluctuation and correlation scales. Indeed, once 
we have fixed the values of any two of the parameters n, x, t on the large scale, the value of 



the third one is uniquely fixed to the same order by the deterministic relation (1.5). Then, 



the random fluctuations of any of these quantities characterize the degree of violation of 



this relation. In other words, we fix a point on the 2D surface defined by the relation 



(1.5) in 3D space of parameters "f,u,v. Then, the small fluctuations in the vicinity of 
this point are represented by an infinitesimal vector normal to the surface, which can be 
projected to one of three directions 7, u, v or any other direction in 3D space. A choice of 
the direction affects only the angle-dependent constants defining the fluctuation scale, while 
the functional form of the distributions is universal. Furthermore, the correlations between 
fluctuations associated with different points of the surface are also universal, as far as the 
points are separated by a distance of order of correlation scale, N 13 . The universality holds 
as the mutual positions of the points vary in a wide range. Indeed, the limiting correlation 
functions of both positions [T6HT8] and times [19] chosen within the domains ( I.ll|1.12 ) and 



( I.13[1.15 ), respectively, yield Airyi correlations for the case of step initial conditions. 

How rigid the universality with respect to the choice of points within the correlation 
function was clarified by Ferrari in [20], whose arguments were based on the observed slow 
decorrelation phenomena. He explained that the limiting correlations can be of two types 
depending on whether the point configurations under consideration are space-like or time- 
like. The correlations for the space-like configurations are, up to a non- universal scaling 
factor, of the same form as the purely spacial correlations. Specifically, when the distance 
between points is of order N°, the fluctuations at these points are described by the Airyi, 
Airy 2 e.t.c. ensembles, depending on the initial conditions, like in the purely spacial case. 
However, if the point configuration is time-like, the fluctuations, measured at the character- 
istic fluctuation scale N a , remain fully correlated, i.e. identical, until the distance between 
the points will be of order of N, which is much larger than N@. 

The definitions of space-like and time-like point configurations used in [20] for the polynu- 
clear growth (PNG) model and extended by Corwin, Ferrari and Peche (CFP), [21], to a 
wide range of other models including TASEP were, however, different from the one accepted 
in [T71 ITS] . To classify our results correctly, we recap here the main idea of CFP. Their 
formulation used the language of the last passage percolation [3], which can be directly, 
mapped to the TASEP as well as to many different models [2T] . Let be the first quad- 
rant of R 2 . Each point of IR+ with positive integer coordinates G N 2 f] R+, is assigned 
a geometrically distributed random variable Tij, 

P(T h3 =t)=q\l-q). (1.16) 



6 



A particular realization of the TASEP evolution is recorded in the values of Tij. Namely, 
Tij is the time the i-th particle is waiting for before making j-th step after it has been 
allowed to move. A directed lattice paths, n^^)-^^), is the path, which starts at the 
point (xi,yi) and, making only unit steps either upward, — > + 1), or rightward, 
— > (i + 1, j), ends at the point (a^jfe)- The sum of Tjj over the path is referred to as 
the last passage time. As it was shown by Johansson for the TASEP with parallel update 
[3], the last passage time, maximized over the set of all paths from (1, 1) to (n, m), 

T n ,m = max V T K j (1.17) 

{n (0 .0)-,(n,m)} (ij . )en 

is related to time t n (m) the n-th particle takes to make m steps, t n (m) = T n ,m + n. For the 
TASEP with backward sequential update these two times are simply equal, t n (m) = T n ,m- 
Other models can be obtained as limiting cases. In the limit q — > 1 with rescaling of time t — > 
t(l — q) we obtain the exponential distribution of waiting times, which defines the continuous 
time TASEP. In the opposite limit q — > the first quadrant is filled mainly by zeroes, 
while "one" appears rarely having concentration q. After going to the continuous limit with 
rescaled coordinates (x, y) — > (qx, qy), the distribution of "ones" on the background of zeroes 
becomes the Poisson process in the first quadrant, which in turn can be used to define the 



PNG [T5l [22] . Given (n, m), the probability distribution of waiting times (1.16) induces the 
distribution P(7^, m < a) of the last passage time 7^ m . The joint distributions P{J~m,m\ < 
ai, . . . , Tn k , mk < Ofc) of the last passage times for k different points (ni, mj), . . . , (n k , m k ) are 
referred to as fc-point correlation functions. 

According to CFP, two-point configuration ((m,mi), (^2,^2)) is time-like if the points 
can be connected by a directed path II( niimi )_>.( n2 >ma ) and is space-like otherwise. Suppose 
that rii < ri2- Obviously, the time-like conditions are 

mi < m-i when n\ < ri2, (1-18) 
m% < m 2 when rii = n 2 . 

Recall that in the TASEP with step initial conditions a particle with the number n starts 
at initial position x° = — n + 1. Therefore, the spatial coordinate of the particle, which has 
traveled for the distance m, is x n = m — n + 1. Then, the space-like condition opposite to 



(1.18) can be translated to the one for the space coordinates: 

Xi < Xj when i > j. (1-19) 



This is the condition that the slow decorrelation does not occur, and, correspondingly, the 



universality holds. One can see that the points (1.14) of final configurations within the 
current correlation functions satisfy this condition. Also, due to non-crossing of particle 
trajectories, these conditions hold automatically when the time moments are chosen in the 
domain ( I.ll|1.12 ). Therefore, the point configurations studied in [T6H19] are space-like 
according to CFP classification. However, in the complementary domain, both types of 
the scaling behaviour present. Thus, the division to time-like and space-like configurations 
proposed by CFP is more adequate if one wants to distinguish between different types of 
universal behaviour of correlation functions. By this reason, we keep on their terminology, 



where the space-like configurations in TASEP are defined by the condition (1.19) and time- 
like by the opposite one. The current correlation functions calculated in [19] were just an 
example of space-like correlations beyond the domain studied in [171 [18]. In fact, the earliest 
result on space-like correlations was obtained in [22J, where the universality of the scaling 
limit was shown in context of the PNG model in the whole space-like domain. However the 
microscopic consideration in context of the TASEP was limited to (|TTT||lT2[ in [IMS] and 



to ( I.13|1.15 ) in [19], where the spacial coordinates were fixed by (1.14). 

In this paper we extend the microscopic derivation of the TASEP correlation functions 
to the rest of the space-like domain, what has not been covered by previous analysis. 



C. General overview and the aim of the present work. 

We conclude the introductory part with an informal outline of the recent development of 
the theory of multipoint correlation functions described above and formulation of purposes 
we are going to fulfil below. 

Though the previous results were formulated in terms of distributions of various quan- 
tities, they can be considered in a similar fashion if we look at the TASEP as at the prob- 
ability measure over collections of interacting lattice paths (the space-time trajectories of 
particles), which can go one step down (particle stays) or down-right (particle makes a step) 
in the space-time plane. Then the correlation functions give marginal probabilities of cer- 
tain points or bonds of the underlying lattice to belong to paths corresponding to selected 
particles. Specifically the development can be roughly divided into three stages depicted in 
Fig. (IT]) . At the first stage the points were fixed at the same moment of time, e.g. those 
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(a) 



(b) 



(c) 



FIG. 1: (a) Equal-time correlation functions; (b) The first extension of the space-time domain; 
(c) Current correlation functions. 

encircled in Fig.([]Ji). The basic achievement of this stage, mentioned in subsection 1.1, is 
revealing the structure of determinantal process in the TASEP [6j [7] . 

The second stage described in subsection 1.2 is characterized by an extension of the range 
of point configurations to space-time domain shown in Fig.Qb)). The condition crucial for 
the solution is the possibility to cut off the part of particle trajectory following the selected 
point without affecting the remaining part. In the first case we just stopped at the moment 
of interest and the independence from the future was a trivial consequence of the fact that 
the TASEP is a Markov process. In the second case similar independence follows from 
another Markov property specific for the TASEP dynamics [T71 HH]: the particles in the 
TASEP do not affect an evolution of other particles to the right of them. Therefore, one can 
drop a part of a particle trajectory if there is no points fixed to the left of it at later time, 
see Fig.(JTTb)), so that the time corresponding to the selected points increases weakly from 
left to right. Finally one again arrives at the determinatal process, though more elaborated 
than the one in the first case. 



The third stage, referred to as current correlation functions, is depicted in Fig. (TT 



Here the particle trajectories propagate equal distances in spatial direction and the selected 
points are fixed at different moments of time, which, as seen from the picture, must increase 
weakly from right to left. At the first glance this situation is in contradiction with the 
above "trajectory cutting" ideology. However it is not difficult to convince oneself that 
if we require that the trajectory makes a step forward after the selected point, it has no 
chance to interact with the trajectory that ends one step to the left of it at later time. 
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FIG. 2: The extension considered in the present work. The time moments are not ordered, 
while the space positions are. The steps of first two particles when they leave the zone of possible 
interaction with the following particles are shown by dotted lines. 



Therefore the part of the trajectory after this step can be dropped. This is a Markov 
property analogous to the previous one, which lies behind the solution. Technically, the 
reduction of the number of particles continuing evolution can be performed by use of so 
called generalized Green functions introduced in [23J and applied in [19], which in turn can 
be reduced to the determinantal process again. On the language of lattice paths this solution 
yeilds the probability of having a fixed bonds within the trajectories selected particles. 

Our goal here is to unify all the previous achievements. Below we calculate the proba- 
bilities of trajectories of selected particles to contain given points or bonds, as shown in in 
Fig.Q. The range of point configurations we consider is wider then in the earlier solutions. 
Combination of two above Markov properties and use of the generalized Green function 
allow us to remove time ordering completely. The tools we use, however, are applicable only 
when the spacial positions of the endpoints are strictly ordered in space. This is the only 



major constraint, which is nothing but the space-like condition described above (1.19) 



Though the ensemble of lattice paths gives a good pictorial representation of the problem, 
this language is not suitable for real calculations and presentation of the results, because 
the whole set of lattice paths is too big. To quantify the results we need a suitable proba- 
bility space, where we could enumerate all our possible random outcomes. In the solutions 
mentioned above this was the set of particle coordinates (Z), i.e. the lower horizontal line 
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in Fig. [jja), product of several such sets, i.e. subsequent horizontal lines in Fig. [jjb), or the 
set of exit times enumerating the points at the vertical lines in Fig. [T|(b) , respectively. Let 
us think about these lines as the boundaries dividing the space-time plane into two parts. 
In all cases the space-time trajectories of particles go from one part to another right at the 
points we select. Therefore we can think of the probabilities under consideration as the 
probabilities for particle trajectories to go from the boundary at specified points. Known as 
exit probabilities such quantities are important in the extremal statistics [21]. Exit proba- 
bilities is a convenient language to represent most general correlation functions. To extend 
the range of space time configurations, we consider the boundaries of more general form: 
a broken line going from northeast to southwest by unit steps either vertical or horizontal, 
which divides the space-time plane into two parts. Consider now the space time trajectory 
of a single particle starting at the northwest part. Obviously, going from the northwest to 
southeast, this trajectory will finally traverse the boundary. The question is, where will it 
happen? We can enumerate the sites of the plane belonging to the boundary by a single 
generalized coordinate r = t — x, which runs over Z. The value of r corresponding to the site 
where the trajectory exits the boundary is a random variable, and its distribution P{r < a) 
is the quantity of interest. The probability distribution of particle coordinate at specified 
time moment and of the time the particle jumps from a specified site are particular cases of 
this general quantity. Note that the exit occurs by two ways (down and down-right) from 
horizontal parts of the boundary and only down-right from vertical parts in the same way 
as above. 

The problem we address below is a direct generalization of one-particle picture described. 
We consider a collection of m arbitrary boundaries, each with its own space-time coordinate 
t% running in Z, and enquire about the joint distribution P(t\ < ai, . . . ,r m < a m ) of the 
coordinates of sites at which specified particles go from given boundaries, see Fig|3j This 
construction allows one to remove any time ordering constraints and include into the scheme 
a possibility to consider both probability of particle being at a site and jumping from it. 
The geometric constraints on the boundaries from which the constraint on the accessible 
point configuration follow will be detailed in the next section. 

After obtaining the results on exit probabilities we perform the scaling analysis of the 
formulas obtained. The lattice boundaries can be used to approximate smooth curves in 
the plane, and the selected points are considered in the vicinity of a smooth path traversing 
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FIG. 3: Exit probabilities. The broken lines are the boundaries. An obligatory step forward must 
be added for exit to happen from vertical parts of the boundaries. The trajectory exiting from the 
horizontal part can continue in any of two ways. 



these curves. The main claim stemming from this analysis is that the large scale behaviour 
of the of exit probabilities is universal as far as the path under consideration do not violate 
the space-like constraint: the fluctuations of generalized exit coordinates of particles starting 
from step initial conditions are described by Airy2 ensemble in the same way as in purely 
spacial case. 

The article is organized as follows. In the section |TT] we give definitions and formulate 
two main results of the paper: exit probability distribution for trajectories of finite number 



of particles at the lattice (Theorem II. 4) and its scaling limit (Theorem II. 5). In section 



III we reformulate the TASEP in terms of signed determinantal process and prove theorem 



II. 4 about exact form of the correlation function. The section IV is devoted to asymptotic 



analysis of the results of previous sections, were we prove theorem II. 5 



II. METHOD AND RESULTS 



A. Exit probabilities for particle trajectories on the space-time lattice 

To define exit probability for a single particle performing ID asymmetric random walk, 
consider a decomposition of the space-time 2D lattice into two complementary subsets 
Q |J Q = 1? . Given the random walk having started at point (a; , £°) £ Q, the exit probabil- 
ity referring to Q is a probability distribution of subsets of the boundary of Q from which 
the particle exits Q. We will consider only sets having a property that once the particle has 
exited Q, it never returns there again. Then the probability of exit from given point of the 
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boundary does not depend on the global form of the boundary of Q. Rather it is simply a 
product of the probability for the particle trajectory to reach this point and the probability 
that the step from this points results an exit from Q. This is the case if the boundary of fl 
is defined in the following way. 

Definition II. 1 The boundary B is an infinite countable subset of 1? 

B = {b(r)eZ 2 } T ^ (III) 

with the following staircase-like structure. Let b{r) = (x,t). Then the next point of the 
boundary will be either 

6(r+l) = (x-l,t) (II.2) 

or 

b(r + l) = (x,t + l), (II.3) 

for any r G 7L. A natural integer variable r increasing along the boundary from north-east 
to southwest can be chosen as r = t — x, (x, t) G B. 

Note that this construction ensures that the trajectory of a particle started in fl eventually 
leaves fl through the points of the boundary B with probability one and never returns there 
again. The probability distribution of the sets of these points is a simplest example of the 
problem we address here. More generally one can consider a collection of embedded sets 
ill C ^2 C . . . , with boundaries Bi, B2, ■ ■ ■ and look for the joint distribution of successive 
exits from these boundaries. 

The idea of exit probabilities for N particles undergoing the TASEP evolution on ID lat- 
tice generalizes the single-particle picture. Now we are interested in how the trajectories of 
collection of interacting particles exit given sets. The quantity of interest is the joint distri- 
bution of subsets of their boundaries at which exits occur. Again, great simplification takes 
place i) for such boundaries, that once the trajectories exited them they never return there 
again. On the other hand we would like that for many particles ii) all possible configurations 
of exit points on the collection of boundaries would be assigned a probability measure in 
the same way as the points of the boundary in single-particle case. The main tool which 
allows us to work with exit probabilities is the Generalized Green Function (GGF). Unlike 
purely spatial Green function used by other authors, the GGF allows us to work directly 



13 



with space-time point configurations (x,t) = ((xi,ti), . . . , (xjv,tjv)) belonging to the set of 
admissible configurations defined by constraints 

X\ > X2 > ■ ■ ■ > Xn (II.4) 
ti < t 2 < • ■ • < t N . (11.5) 

For iV particles the concept of the boundary can be generalized to iV-boundary, which 
allows us meet (i) as well as (ii). 

Definition II. 2 Given boundary B, the N -boundary 3n C {1, . . . , N} x is defined as a 
disjoint union of N copies of B, 

N 

B N =[_\ B k , (II.6) 

fc=i 

where the copy B^ = {bk{i)}iew, 2 associated with k-th particle is shifted by (k — 1) steps back 
with respect to the first one in horizontal (spacial) direction of space-time plane, 

b k (i) = (x(i)-k + l,t(i)), (II.7) 

fc = l,...,JV. 

The iV-boundary is a generalization of the line with fixed time coordinate and of the 
set of lines with fixed space coordinates, which where the probability spaces used in [TTl 
[18] and in [19] respectively. Having started from an admissible point configuration, iV 
particle trajectories will reach given iV-boundary after some evolution, traverse it and go 
from some points of the iV-boundary to continue the evolution. Then, the non-crossing of 
the trajectories ensures that the configuration of the departure points at the iV-boundary is 
admissible as well. 

To specify from which to which point sets the system can pass in course of the TASEP 
evolution, we also need a relation between subsets of {1, . . . , iV} x 1? . 

Definition II. 3 Let ft, O' C {1, . . . , iV} x 1? . We say that relation 

tt^tt' (11.8) 

holds, if for any (xk,tk) G Q and any (x' h ,t' k ) G fl' 

(4, t' k ) e{(x,t):t> t k } \J{(x, t):x> x k }. (II.9) 
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FIG. 4: Trajectories of five TASEP particles traversing three 5-boundaries. Black segments em- 
phasize that particles make a compulsory step forward at the sites belonging to vertical parts 
of boundaries from which the exits occur. The exits included into correlation function with 
Ni = N 2 = 5,N 3 = 4,A 4 = 3,iV 5 = 2,N 6 = 1 and jfej = 1, k 2 = k 3 = 2, k A = k 5 = k 6 = 3 
are shown in circles. 

Note that the subindices denote the variable from the set {1, . . . , N} and are associated 
with the number of a particle. 

As it was explained in [19], a space-time trajectory of a particle starting from a point 
preceding to a given boundary, eventually transverses the boundary with probability one. 
The question we address is: What is the probability for the trajectory to go from a given 
subset of the boundary? More generally we address the same question to a collection of 
particles and a set of points at several boundaries. 

To be specific, consider the TASEP evolution of N particles governed by the dynamical 
rules I-III. Let the initial configuration x° be defined by 

a£ = -i + l, i = l,...,N. (11.10) 

Let us fix a collection of A^-boundaries, B , . . . , B m , m > 0, such that 

x° -< B 1 -< < B m . (11.11) 

and fix the one-particle boundaries , . . . , £>^ ; within the A"-boundaries. Here the upper 
indices 1 = k± < ■ ■ ■ < hi — m refer to the number of A^-boundary, the lower indices, 
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N < ■ ■ ■ < Ni = N, to the particle number, and I > m. We suggest that at least one 
particle is fixed at each iV-boundary, i.e. either k i+ i = h t or k i+ i = hi + 1. We also 
require that equality Ni = N i+ i for some i suggests that k i+ i = hi + 1, i.e. two subsequent 
space-time points chosen for one particle should be put onto subsequent ^-boundaries, 
and no other particles with number less than N can be fixed at the iV-boundary fcj. Let 
space-time positions of points b*(i) within the corresponding boundary B k be indexed by 
index i G 7L in the same way as in Defs. ( U.1[II.2 ). The quantity of interest is the joint 



probability distribution P[i\ < a\, . . . , %i < ai) of the points (6^ (ii), . . . , from which 

the space-time trajectories of particles Ni, . . . , N make steps when leaving the boundaries 
(jBJ^, . . . , B^) respectively. 

The first main result of the present paper can be stated as the following theorem. 

Theorem II. 4 Under the above conditions the joint probability distribution of exit points 
is given by the Fredholm determinant 

P(ii < ai, . . . , %i < a t ) = det(l - r] a Kr] a ) l2{{B ^ ^ }) (11.12) 



with the kernel 



K(b%,b%)= (11.13) 



dv f dw (l-p(^)P (v-l) N i v Xn j 



1 2tcw J To v 2mw (w — v) (1/v + l/7r 2 — 1) 



1(N 2 > Nx 



T To i 2mw 2 (w - l) N s~ N *(l/v + 1/tt 2 - 1)' 

where r] a = > a x ) x ■ ■ ■ x l(ix > a m ), b k j. = (x Ni ,t Ni ) G B k £., i, j = 1, . . . ,1 and tt 2 = l,p 
is the probability of step from the boundary B^. at point b^.. 



B. Scaling limit of correlation functions 

In the large scale the boundaries can be treated as approximations of continuous differ- 
entiable paths in the space-time plane. Consider a scaling limit associated with sending 
to infinity a large parameter L — > oo, as the time-space coordinates and particle numbers 
measured at L-scale are fixed: x/L, n/L, t/L = const. Let us introduce variable change 
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r = t — x 
t + x 



Lx 



(11.14) 
(11.15) 



As it was noted earlier the variable (11.14) naturally enumerates points at the boundary. 



Correspondingly, the function ((x, 9) defines a one-parameter family of curves spanning the 
whole space-time plane as 9 varies in R. As the parameter x runs in R, it defines a point at 
a particular curve corresponding to some fixed value of 9. The properties of ((x, 9) follow 
from the properties of boundaries. Specifically, we suggest that 



and 

We now suppose that for k = 1, 
sponding to fixed set (#i, . . . , 9 m )\ 

b\{[Lx]) = L. 



d_ 
89 



dx 
d_ 



< 1 



C(X,0)>1. 



(11.16) 



(11.17) 



. , m the boundaries B\ approximate the curves corre- 

C(xA)-x C{x,o k ) + x 



+ o(L°) 



(11.18) 



2 ' 2 

where the notation [ ] is for integer part of a real number and the correction term should 
not contribute on a characteristic fluctuation scale, i.e. a = 1/3. For technical purposes 
we will suggest that the correction term is uniform over the boundary. These boundaries 
correspond to the first particle. For general particle with number n = [Lu] we have to 
consider the boundary B\ shifting the spacial coordinate by n — 1 steps backward: 



b k n([L X }) = L 



(11.19) 



Recall that on the large scale, > oo, the trajectories of particles are 



deterministic, defined by the relation (1.5). In terms of new variables the relation turns into 

Vp(C(x, W+xj- V(C(x, o)-x)- = 0, (11.20) 

which uniquely fixes value of x given those of 9 and u, provided that the corresponding curve 
passes through the rarefaction fan defined by 



X<C(X,9)<\ ±Z X- 
1 — p 



(11.21) 
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Let us consider a path in 9 — v plane: 

9 = 9{r) 
v = u(r) 

with different iable functions 9(r) and v(r), such that 



r e R, 



(11.22) 



89 

dr 



>0, 



dv 
dr 



< 



(11.23) 



and 



dv 



> 1. 



(11.24) 



09 

dr dr 

We select m points at the path, r = r 1; . . . , r m , so that the integers Ni, . . . , N m from Theo- 



rem 1.1 are given by iVj = [Lz/(rj)], and 9i = 9(ri). The inequalities (11.23) and (11.24) then 
guarantee that the constraints on k\, . . . , k m and N\, . . . , N m from Theorem 1.1 are satis- 
fied and together with non-crossing of particle trajectories ensure that points of this path 
accessible for particle trajectories with nonzero probability form space-like configurations. 



Substituting functions 9(r) and u(r) into (11.20) we obtain an equation, which, given 
r, can be resolved with respect to \. For a given path a unique solution exists for any r 
within the range, in which the boundary corresponding to 9(r) passes trough the rarefaction 



fan (11.21). This solution is a monotonous function of r, which we denote x( r )- It defines 



the macroscopic deterministic location of the point, from where given particle exits given 
boundary, see Fig. [5| We are now turn to the fluctuations of these points referred to the 
boundaries and particle numbers separated by the distances of order of correlation length 
from each other. Suppose that 

n = r + Ui L- 1/3 . (11.25) 

The corresponding values of \ are given by their deterministic parts x( r i) P ms a random 
variable of order of fluctuation scale 



Xi = X{r 



) ■ 0. (»-!/. 



(11.26) 



In what follows we show that the random variable converges to the universal Airy 2 



process for a class boundaries, which can be approximated by (11.18 )-(!!. 19). 



Theorem II. 5 The following limit holds in a sense of finite- dimensional distributions: 



lim £ L = KfA 2 (K c u), 

L— >-oo 



(11.27) 
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N l N 2 N 3 




FIG. 5: Exit probabilities on a space-like path in x-t plane. The wedge bounded by black straight 
lines is the rarefaction fan area. The deterministic trajectories of particles with numbers N\, N2,N^ 
are shown in red. The green lines are boundaries with coordinates x = L((£(#j,x) — x)/2 — Vi) 
and t = L(Q(6i, %) + x)/2, where i = 1, 2, 3, corresponding to three fixed values of 9: Q\ < 62 < #3. 
Dashed line is the projection of the path {v(r),Q(r)) to x-t plane: x = L((((r) — x(r))/2 — v(r)), 
t = L(C(r) + x( r ))/2- The black dots are the points where exits occur. 

where A2 is the Airy 2 process characterized by multipoint distributions: 

PTob(A 2 {u 1 ) < s u . . .,A 2 {u m ) < s m ) 

= det (1 - VsKMry 2 Vs)) L 2 ({ni> ... jnm}xm • (H.28) 

where in the r.h.s. we have the extended Airy kernel, 

K A ny 2 (6,Ci;aC 2 ) (11-29) 




/ °° d\e x ^Ai(\ + Ci)Ai(A + Ca), 6 < 6 
-/_° oo rfAe A ^-^)Ai(A + Ci)Ai(A + C 2 ), 6 > £1 ' 



The model dependent constants k c and Kf defining the correlation and fluctuation scales 
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respectively are given by 



P 1/6 (v^-v^7)- 1/3 (v^-v^)' 1/3 m 
c 271/60,1/6 ( v ^(i + C (o,i)( ro )) + v ^(i_ C (o,i)( ro ))) 1 • > 



qv' (r ) (C (r ) - X (n>) C (0,1) (r )) (y/pu ~ jl)' 1 
- 9' (r ) C (1 ' 0) (r )(^w-V7)] 



(^*J - V7) 1/3 (VP7(1 + C m (r )) + V^(l - C (0 ' 1} (ro))) 
K / = ^- ^273 ' ( IL31 ) 



2p 1 / 6 W 1 / 3 7 1 /3 (^ - ^y) 

where we denote £(r) = £(0(r), x(r)) ; C ( ' 1 ' ' ) ( r o) r*C < ~°' 1 ' ) ( r o) derivative of the function 
((9,x) with respect to the first (second) argument at the point (9(tq), x( r o)) and parameters 
7 and u are those defined in (1.2), 7 = (C( r o) ~ x( r o))/2 and u = (((ro) + x( r o))/2. 

The non-universal constants Kf and k c are the most general ones for the TASEP with 
backward update. They depend not only on the macroscopic space-time location defined by 
C(r ) and x( r o)> but also on the local slope and local density of the boundaries at this point 
via the derivatives C ( ' ' 1 ' ) ( r o) an d C < ~ 1 '°' ) ( r o)) respectively. Particular cases studied before can 
easily be restored from the expressions obtained. For example, for purely spacial boundary 
used for measuring particle coordinates at fixed time we can take ((6, x) = 2t — x, while the 
case of current correlation functions [19] corresponds to ((6, x) — 2x + %■ F° r the space-like 
correlation functions of particle coordinates studied in [T71 [18] we take ((9, x) = 29 — x, and 
the tagged particle case [16] corresponds to i/'(r) = 0. 



III. DETERMINANTAL POINT PROCESSES ON THE BOUNDARIES 



A. Single iV-boundary. 

We first introduce the Generalized Green Function (GGF) using the determinantal for- 
mula proposed in [23] and proved in [19| . which generalizes the formulae of simple Green 
function obtained in [26] for continuous time TASEP and generalized to the backward se- 
quential update in [25J 

Given two admissible configurations 

b° = (x ,t ) = ((xU°i),- 
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and 

b = (x,t) = ((x 1 ,t 1 ), . . . , (x N ,t N )) 

we define 

G(b\b°) = det [Fj_i(bi - 6°)] ^ N . (III.l) 
where {pi — 6°) = (xj — x®, ti — t®) is componentwise extraction and 

f ^ <L — (q+ R Y (1 - w)- n w x , t > 
F n (x,t)=^ 2 * ljTo w V ™ ; V ; . (III.2) 

jO, t<0 
For point = (x, i) e jB at the boundary, we introduce an exit probability 

I p, if Hi + l) = (x,t+ 1) 
vr 6 (6(0) = \ P ' 1 J 1 \ (III.3) 

[ 1, if 6(i + l) = (x + l,t) 
and for iV-point configuration b E B 



N 

n B (b) = l[n^(b k ), (III.4) 

k=i 

where the subscript k specifies a boundary within the iV-bounday, or the associated particle. 
The function 

g(b\b°) = 7i B (b)G(b\b°) (ins) 

gives probability for the space-time trajectories of particles to go away from the boundary 
via the points of b, given they started from b°. 

We now show that this probability can be reinterpreted in terms of an auxiliary signed 
determinantal point process on B. Consider a signed measure on Z> T() x {1, . . . , iV}, 

JV-l 

M en = z~ II det ^> Tf 1 )]^ det [K-At? )]y =1 , (me) 

W n=0 

assigned to the sets of the form 

r= U {C<Ci.<-,<flc% x{i,...,4 (IH.7) 

l<n<iV 

Here we define the function 

, Tr Bn+1 (b n+1 (y)), y> z 
M*,V) = { ^W'V- (lug) 

0, y < z 
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for x, y E Z> T0 and bk(y) E Bk, k = 1, . . . , N and the function 

9Z{t) = (-l) h F- k (b N (t)-b%_ h ), (IIL9) 

where 

F n (x,t) = ^{ — fg+ -V(l - ^)—^. (111.10) 



27ri 7r w ^ ^ 



The integral representation holds for t E 7L. This is unlike F n (x,t) which coincides with 



F n (x,t) when i > and vanishes at £ < 0, see (III. 2). The numbers r/ are integers bounded 
by number To from below. The number To is chosen so that (tq) = 0, which is always 
possible by construction of the boundaries. 

We also introduce fictitious variables t™ -1 , 1 < n < N, which are fixed to t™ _1 = tq. 
Thus for any r- E Zj> to we have 

0n(r n " +1 ,r™ +1 ) ee vr^(W^ +1 )) ( IIL11 ) 

for j = 1, ..,n + 1. The numbers r™, i = l,...,n are mapped to the sites b n (r™) on B n . 
Therefore the measure on the iV-boundary 13 is naturally defined as pushforward of A4(T) 
under this mapping. 

One can consider 6 n (r"), 1 < j < n < N, as coordinates of auxiliary fictitious particles 
indexed by j leaving at the boundaries B n . These particles evolve as shown in FigJSJ First, 
N particles arrive from their initial state encoded in the functions (t) at the points of 
boundary Bn with numbers rf, . . . , rjy > r . Then, they jump to the sites of the boundary 
Bn-i with the same numbers, and go up along the boundary £>at_i (from south-west to 
north-east, so that the number r indexing position at the boundary decreases) any distance 
respecting mutual noncrossing of particle trajectories. The weight of the jump between the 
boundaries, outgoing from a site b, is ir BN (b) and the weight of going along the boundary 
is 1 independently of the distance. The last (iV-th) particle is forced to go to the resevoir 
(r^ -1 = Tq) and disappear. The final positions of the other particles at the boundary B^-i 
are denoted t^ 1 , ■ ■ ■ , t^Zx, from which they jump to the boundary B^-2 with the weights 
tt Bn - 1 (ti~ 1 ), . . . ,k Bn ~ 1 (tnZx)i e.t.c. The process is repeated until the particle number 1 
jumps from the point 6i(r 1 1 ) of B\ and disappears. This picture generalizes the auxiliary 
processes described for the cases of constant time [6j [H [T7J [18] and fixed spacial coordinates 
[19J to the case of general boundaries. 
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<l>2 ^1 



FIG. 6: Evolution of fictitious particles for N = 3. The first stage from initial points to the 
boundary B% is described by the functions ( r & )• The following stages from B% to £>2, from £>2 
to B\ and from B\ to Bo are encoded in ^2 > </>i and ^>o respectively. The coordinates of real TASEP 
particles are shown by black dots. 



The fictitious particles are similar to vicious walkers (or free fermions), which can be 
seen from the Karlin-McGregor-Lindstrom-Gessel-Viennot [27H29] determinantal form of the 



transition weights entering the product (III. 6) that ensure nonintersecting of their space 



time trajectories. The last determinant can be treated as integrated with given initial 
distribution. Such a free fermionic structure allows calculation of the correlation functions 
for fictitious particles, which turns out to be determinantal. On the other hand, below we 
show that the joint distribution of N positions of the first fictitious particle obtained by 



integration of the measure (III. 6) over the positions of the other particles coincides with the 



Green function of TASEP. Thus the problem of correlations in TASEP can be reduced to 
calculation of correlations between noninteracting mutually avoiding fictitious particles. 
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To show that the GGF can be interpreted in terms of the measure M. (T) we prove the 
following proposition. 

Proposition III.l Given N -boundary 3, initial and final configurations b° -< 3 and b C 3 

respectively, GGF Q B (b\b°) associated with the boundary 3 is a marginal of the measure M.: 

g B (b\b°) = m Qj{7f = i k yj , (HI.12) 

where determine the location of the points of b = (&i (»i) , ■ • • , &jv(*jv)) a t corre- 

sponding boundaries within the N -boundary. 

To prove this statement, one represents the GGF as a sum over the boundary points 
in a way similar to that used for space variables in [6], [7J H7J [18] and for time variables in 
[5j [19]. The proof of the summation uses contiguous relations for the values of the function 
F n (b(r)) at adjacent points of the boundary, which unify similar relations for space and time 
variables. 

Lemma III. 2 Let b k {r) be the point at the boundary B k within the A? -boundary 3. Then 
the contiguous relations hold for the function F n {bk(r)) 

7T B (b k (T))F n (b k (T)) = F n+1 (b k ^(r + 1)) - F n+1 (b k ^(r)) (111.13) 

Proof. The relation to be proved is in fact two contiguous relations for the function F n (x, t), 

F n (x,t) = F n+1 (x,t)- F n+1 (x + l,t), (111.14) 
P F n (x, t) = F n+l (x + 1, t + 1) - F n+1 (x + 1, t), (111.15) 

as one relation. The two latter relations follow from the integral representation of the 
function F n (x, t). m 
Then we have: 

Lemma III. 3 Given N -boundary 3, initial and final configurations b° ~< 3 and b C 3 

respectively, the function Q(b\b°) can be represented as a sum: 

$( b \ b °)) = e n *** (w)) 

D l<i<n<N 

x (-1)^ det^^MT^) - &Sr-i)]JB> ( IIL1 6) 
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where bj(r[) = bj, j = 1, . . . , N, and the summation variables take their values in the domain 



D = {rle Z> T0 ,2 < i < j < N\ij > rf\rl > r^} 



(111.17) 



Proof. Using the contiguous relation (III. 13) the proof just follows the similar proofs in 



[HHTJ IT7HTU] . Note that the lower summation bound To is chosen such that the functions F 
under the determinant vanish at this point. Indeed this is true for F n (x,t) when x > t. By 
construction of boundaries it is always possible to find suitable r to ensure this inequality. 
To be specific we choose the maximal of these numbers. ■ 



To complete the proof of proposition ( III. 1 ) we need to show that the summation over 



the domain D can be replaced by the summation over the sets of the form (III. 7) 



Lemma III. 4 The domain of summation in (III. 16) can be replaced by 



D = {rle Z> T0) 2<i<3<N\Tj<4_ x } 



(HI.18) 



Proof. Apparently inequalities in (III. 17) suggest those in (III. 18). We also need to show 



the converse: the measure (III. 6) is zero everywhere in D unless the inequalities from in 



(III. 17) are satisfied. The statement can be proved by reproducing the arguments from [19 



To find the correlation functions of the TASEP we first calculate the correlation functions 
of the measure M.(T). The functional form of Ai(T) suggests that the correlation functions 
are determinantal. Derivation of the correlation kernel was explained in great detail in [7J. 
To proceed with the calculation, we introduce convolution 



<P {ni ' n2) (x,y) 



h * 0m +i * ' ' ' * 4>n 2 ~i)(x,y), n x <n 2 



(111.19) 



where (a * b)(x, y) = E* e z> TQ a ( x i z )K z i v)i and 



^»*j^)(r). 



(111.20) 



Note that in terms of the coordinates of fictitious particles function (fi( ni ' n2 \x,y) is 
the transition weight between points at the boundaries B ni and B n2 . Hence, the points 



parameterized by the variables x and y in (III. 19) live at B m and B n2 , respectively, while 



the argument of ^-j( T ) ^ n (HI. 20) lives on the boundary B n . 
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Consider functions 



{(0o * 4> {l ' n) ){rlrl . . . , (0 n _ 2 * r^lCi.^^il^.T)}. (IIL21) 

They are linearly independent and hence can serve as a basis of an n-dimensional linear 
space V n . We construct another basis of V n , {$™(r), j — 0, . . . , n — 1}, which is fixed by the 
orthogonality relations 

53$?(r)«y(r) = tfy. (HI.22) 

Then, under the 

Assumption (A) : 0n(r™ +1 , r) = c„$q (r) with some Cn. ^ 0, n = 1, . . . , N, 
the kernel has the form 

ri2 

K(n u r l5 n 2 , r 2 ) = -^'^(n, r 2 ) + £ ^_ fe (Tx)^_ fc (T 2 ). (111.23) 

fc=i 

Applying repeatedly the convolution with 4>n -ii ■ ■ ■ ,<fij to ^^_j(r) we obtain 

Lemma III. 5 Given N -boundary B, the functions \P"(t) have the following integral repre- 
sentation. 

*2(r) = jf ^ -p^^y n(T) - l)*^^)-^-*" 1 , (111.24) 

TTie contour of integration Tq^ encircles the poles w = 0,1, leaving all the other singularities 
outside. 

To find basis V n , we have to specify initial conditions. For usual step initial conditions, 
the orthogonalization can easily be performed. 

Lemma III. 6 Given step initial conditions, b° k = (— k+1, 0) for k = 1, . . . , N, the functions 



V&fcCr) and $™(t) satisfying (III. 22) are given by 



p 1 27ri \ w J \ w 



$»(r) = 4 — ] ^ J! (IH.26) 

' 2vri V v J ((l/ir B (b(r)) - l)u + 1) ' V ; 



where the contour of integration T 1 encircles the pole v = 1 anticlockwise. 
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Proof. The function ^(t) is obtained from (III. 24) by an explicit substitution of the 
step initial conditions. To prove the orthogonality conditions (III. 22) one must evaluate 
the sum ^ TgZ \?£(t)<&!?(t). This is done by an interchange of summation and integration. 
After successive summing the geometric progressions for space-like and time-like parts of 
the boundary and taking into account the pole at v — w, we obtain the desirable result. To 
provide the convergence of the resulting sum we note that the choice of contours ensures 
convergence of the sum for r — > oo, while at the lower limit the sum is truncated at r = s, 
so that x n (s) —t n (s) — 1 = 0. Obviously ^(t) = for r < s because no poles remain inside 
the integration contour for k > 0. ■ 

Note that the form of $™(t) depends on whether the site b(r) belongs to time-like or space- 
like part of the boundary, which is reflected in the term containing the exit probability in 
the denominator. Now we note that the assumption A is fulfilled, 

*S(t) = n B "(b n (T)) = Mr:+x,r), (IIL27) 



and we can write the kernel. The summation in (III. 23) yields 



T,KU(n)KU(T2) (ni.28) 

k=l 

(l-p( 3 ^)) t "i (Tl) w *m(ri) 

dv f dlU (l-p( 1 ^ i ))* n 2( T 2) (f-l)"2 ^n 2 (n) 



l 2ti\v J rQv 2Triw (w — v)(l/v + l/n 2 — 1) 



where 7rf = tt b (6(t2)). 

Observe that the function <f) n (x,y) can be written in the form 

, & (l- p (ti))*»(nK + i(^)^(n)-, t , +1 N 

Mn,r 2 )=(p f 7T . 111.29 

/ r0il 2TT1W 2 [W - 1)(1/V + 1/7T2 - 1) 

After a few convolutions we have 

(ni,n2) (n,r 2 ) = i(n 2 >m) (IIL30) 

r dlV (1 — p(=^)) t ' 1 i( T1 )-*n 2 (r 2 ) u; x ni (ri)-a: ni (T2) 
X J rQ1 2mW 2 (W -l) n 2- n l(l/v + 1/7T2-1) 

Then we obtain: 
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Proposition III. 7 The correlation kernel of the measure M., (III. 6), is 



K(ni,T 1 ;n 2 ,T 2 ) = 

(l-p(^)) f " l(Tl) (1D-I)"l ugniCnj 

(IV f dw (l-p(Hnl))*n2( T 2) {v~l) n 2 ^n 2 (ri) 



(111.31) 



Ti 2niv J Tov 2ti\w (w — v)(l/v + l/n 2 — 1) 

d w _ p(?±ti))*n 1 (ri)-in 2 (-r2) u ,x ni (ri)-2.- ni (r 2 ) 



l(n 2 > ni) 



ro,i 



2iriw 2 



(w - l) n 2- n i(l/w + 1/tt 2 - 1) 



where tt 2 = tt s (&(t2)) and x n (r) = x(t) + n — 1. 



Determinants of the above correlation kernel yield the correlation functions of the measure 



A4, i.e. probabilities of point sets T, (III. 7), having any given subsets. Then, using the 



inclusion-exclusion principle, we can write down joint distribution 

P = M (r D K 1 < ax} Pi ■ ■ ■ f|K m < «»}) , (111.32) 

of sequences {rf 1 , . . . , t^} for any fixed collection 1 < ni < • • • < n m < N, where 1 < m < 
N, in the form 



P = det (1 - XaK Xa )) lH{ni _ nm} 



xZ) 



(111.33) 



where Fredholm determinant is defined as a sum 
( 1 \ n m 

P = E^ E E • • • E det{if k 

n>0 ii,...,i„=l ri>a il r n >a in 

and Xa( n «)(^) = l(i > a «)- This distribution is the TASEP correlation function of interest, 



T k] n ip T j)}kj=l 



(111.34) 



P = Prob ({r ni < ai} p|{r„ 2 < a 2 } p| • • • f]{r nm < a m } 



(111.35) 



and (III. 33) is a particular case of the Theorem II. 4 applied to the case of single iV-boundary. 



Remarkably, the GGF allowed us to treat very wide range of space-time point configurations 
"in one go" , in the same way as the fixed time and space cases were treated in [7] and (TH] , 
respectively. Any admissible point configuration can be processed in this way, when put to 
a suitable boundary. The set of admissible configurations, however, does not exhaust all the 



possibilities. It turns out that the time ordering constraint (II. 5) can also be removed. To 



this end we apply a multicascade procedure, similar to that used in [T7], to a sequence of 
iV-boundaries. 
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B. Multiple iV-boundary case. 



Consider m mutually distinct Af-boundaries, 13 1 -< • • • -< 3 m . In this section we derive the 
joint Z-point probability distribution for positions b^ at which the trajectories of particles 
Ni depart from boundaries B^., where i = 1,2, and the indices Aj and ki satisfy 



the assumptions of Theorem II.4 We suggest that the space-time points within each N- 
boundary are indexed independently by the indices Tj(fc), where the subindex % = 1, . . . , N 
stands for the the number of the boundary £>f within the iV-boundary TS k and the argument 
k = 1, . . . ,m indexes iV-boundaries. According to Defs. II. 1[ II. 2 we first independently 
define an indexing order b\(ri(k)) for the first boundaries within each A^-boundary and then 
translate it to other N — 1 boundaries by the corresponding left shifts. 

Given a fixed collection of integers a%, . . . , ai, we are looking for the joint probability 

P = Prob (fl{bft ( = b%(r N XkM=i ■■ {tnM < a^UbA (111.36) 

for trajectories of particles Ni,...,Ni to leave corresponding boundaries via points b^., 
located above (in terms of the corresponding indices Ti(k)) the sites b 1 ^ (ai). 

Similarly to the case of single A^-boundary, our strategy is to represent this probability 
distribution as a marginal of a signed determinantal measure on a larger set. Suppose that 
the set (hi, . . . , hi) is of the form 

k Pl = ■■■ = k V2 _ x = 1, (III. 37) 

k P2 = = k p3 -i = 2, 



^Pm — • • • — ki — m, 

which defines a collection of m integers 1 = p\ < • • • < p m < I. The quantity of interest can 
be given in terms of measure A4(-) on point sets 

r= U r(k), (HI.38) 

l<fc<m 

where 

T(k)= □ X n (k) (111.39) 

N Pk <n<N Pk+1 



29 



and 

Xn(k) = {r:(k), < r^k), <...,< rr(fc)} C Z> ro(fe) . (HI.40) 

Then for 1 < k < m the collections of integers (r™(k), . . . , r"(fc)) define point configurations 
(fe"(r"(A;)), . . . , &i(t"(A;))) C <B T fc t , which can be treated as coordinates of fictitious particles 
similarly to the single A-boundary case. The pushforward of the measure Ad under this 
mapping is a measure on the collection of A-boundaries B 1 , . • • , B m . An explicit form of 
this measure is as follows. 

M(T) = A/--Met[<^(rf (l))]i< w <^ (IIL41) 
x n[det[^-i(r^(i),r^(i - l))]i< M <^ 

JV P . 

n=7V Pi-i +1 

where we define functions 

^(r^i^r^ii - i)) = *4r\bir\T?«(i - 1))) 

xF (b\(r^(i))-br(r^(i-l))l 

^HK + i(r]; +1 m, ^(0 > r/»(<), 

0, r£ +1 (i) < 7f (<), 

A/" is a normalization constant and for 1 < i < m 

bl'\r k (i - 1)) = (x k (i - 1) + ^ ^U(i - 1) + (HI.43) 

1 — p 1 — p 

Lower cutoff r (k) is separately chosen for every A-boundary B& in such a way, that any 
transitions to these points have zero measure. Specifically, as in the single A-boundary 
case considered in the previous subsection ^^(r) = for any r < r . In addition 
J r i i j_i(r(i), To(i — 1)) = for any r(i) < r (i). Correspondingly, the auxiliary variables 
T^-\k) are fixed to T^~\k) = r (k). 

The relation between the correlation functions in TASEP and the measure Ad is given 
by the following proposition. 



(111.42) 



(0,^(0) 
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Proposition III. 8 Consider the TASEP evolution starting with the initial conditions b° = 
((xi,t1), . . . , (x° N ,t° N )), where b° is admissible configuration. Consider also m mutually dis- 
tinct N -boundaries , b° -< B 1 -< ■ ■ • -< B m . Let (JVi, . . . , iVj) and (&i, . . . ,ki) be collections of 



integers satisfying assumptions of the Theorem II. 4 Then, the joint probability for space- 
time trajectories of particles Ni, . . . , Ni to go from N -boundaries S fcl , . . . , B kl via points 
, . . . , , respectively, given the trajectories of all particles started from the point config- 
uration b°, is a marginal of the measure Al(7~) of the form 

P (\Ji h % = =m(td \J{T^(h) = n(k), l < i < l}j . (111.44) 

Proof. We first note that instead of the iV-boundaries B , . . . , B m we can consider auxiliary 
N p . -boundaries B l Np = (B{, . . . , B l Np ) for i — 1, . . . , m. This is possible because in TASEP 
the trajectories of particles Pi+1, ■ ■ ■ , N do not influence the trajectories 1, . . . ,Pi+l, and no 



point of the former group is fixed after (and on) B % within the correlation function (111.44). 
Given trajectories 1, . . . ,pi, the sum over all realizations of the trajectories pi + 1, . . . , N 
amounts to one. Therefore, after B l has been passed we can drop the former evolution and 
consider only the latter. Thus, we first consider the transition of N\ particles from 6° to B l Np ^ , 
then the transition of N P2 particles from B Np to B% p , e.t.c. (see Fig 4). The probability 
of each transition is given by corresponding iVp.-particle Green function. To ensure the 
admissibility of particle configurations within the Green function and keep its probabilistic 
meaning we require that after each transition the particles do leave the boundaries. This 
suggests that we insert a compulsory step forward at the points belonging to vertical parts 
of the boundaries. To this end, we supply each step of this kind by the factor of p and define 



the starting points for every transition to be of the form (111.43). Finally, the probability of 
interest, ^({^(^(/Ci))}-^^ ), is the following: 

P (U=i{b^ = b%(n(k))}\b°) 
= E fl G ({H{n{j))}ii\{bV{<j - !))}£) x i^W 



(111.45) 



A j=l 

i 0,1 



where {b^ (rj(— 1))}^Li = b° and the summation is over domain 

A = {rj(i) E 'L>T Xj (i ) ,Tj{i) > Tj-i(i), l<j< N Pi , l<i<m} 

\{r Ni (ki)A<i<l}. 
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Using the determinantal formula of the GGF ( III. 1 ) , we have 
P({b%(r Ni (h))}U\b°) 



E 11 7r B > fe )det Fi-j mrjim-bfinik)) 

A k=l L v 



(111.46) 



J l<i,j<N Pk 

In what follows we are going to introduce auxiliary variables r/(fc) in the same way as we 
did for the case of single iV-boundary, with the only difference that there is a separate set for 
every iV-boundary, indexed by an extra argument k. To proceed further we define several 
domains of summation in these variables: 

A = {7^(0 G Z> T0(J) , 1 < k < j < N Pi \ri(i) > ri +1 (i)} (111.47) 

A = E Z> ToW , l<k<j< N n+1 - l\ri{%) > ri +1 (i)} (111.48) 

D* = A \ WW, N Pi < j < N Pi+1 }; D* = D* \ A (HL49) 

D = UT =1 D* (111.50) 

where we set p m +i = 1 + 1 and r*(z) = Tfe(z) for A; = 1, ... , i\T Pi and z = 1, . . . , m. 

Now we apply Lemmas III. 3 and III.4| to each determinant under the product in r.h.s. of 



(III. 46) to represent it as a sum over the auxiliary variables: 



det [F^bKnii)) - b°/(r k (i - l)))]i< w <^ 



det 



'-N Pi+ k(b%Jr^(i))-b°/(r^i-l)) 



D* 



l<k,l<N Pi 



(111.51) 



n=0 



The endpoints, b^{r), of part of the trajectories within a transition between two N Pi 



boundaries are related to the starting points, b®' +1 (t), of the trajectories within the next 



transition by (111.43). The sums over the range of these positions can be evaluated along 



with a few sums in auxiliary variables coupled to them (see Fig|7]): 



£ det[F_ Npi+k (b% p Xr^^) - b^(i - 1)))] 



Kk,KN v 



Di 



x [I det[0 n (7£(z-l),r; 

n=0 



n+li 



1))] 



l<fc,«<n+l 



(111.52) 



det[F (b\(r^ (0) - ^(r^(i - 1)))]i<,,k^ 
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Xc Xa X-j X 








x?(2) 



FIG. 7: Example of summation over the auxiliary variables t^(z) for two ^-boundaries in the case 
N\ = 5 and N P2 = N 2 = 3. The dotted lines show the variables, over which the summation can be 
explicitly performed. The points shown by circles are fixed within the correlation function. 

The last identity can be proved by repeatedly applying formula 

^7T B (6 fc «)F n (c -6 fe (z)) (111.53) 

: F n+1 (c - b k -i(i 2 + 1)) - F n+1 (c - b k -i(ii)), 



which is another form of Lemma III. 2 where Cq = (xo,to) is a pair of arbitrary constants. 
The resulting expression for the joint distribution is 

P({b%(r Ni (h))}Ub°) 

uv H r_ Nl+l yu Ni yi k yL)j-t 0] 

D 



X 



X 



const det[F_ JVl+ ;(6]y 1 (Tf 1 (1)) - b G l}1 < k j< Nl 
ndet[F (&i(r^(z)) - bf{r^{i - l)))]i< M <^ 

i=2 
N Pi -l 

n det^r^),^)^ 



(111.54) 



Jl<fc,Kn+l, 



n=N, 



Pi+i 



33 



where we put N Pm+1 = 0, D m = and therefore D* m = D* m . m 

If we again appeal to the correspondence with coordinates of fictitious particles, we see 
that the indices r*(i) define coordinates of particles at the boundary B l n . The functions F () 
under the product in ( |III.54 ) describe the transitions between two subsequent iV-boundaries, 



while the functions 0(r™, r™ +1 ) are responsible for transitions between subsequent n-th and 
(n + l)-th boundaries within the same iV-boundary. Note that after we summed out part of 
coordinates, some boundaries fell out of the consideration and only the following remained: 

) • • • > °N P2 > °N P2 ) • • • > n N„„ > °V B „, > • • • J D l • 



Therefore, it is convenient to develop another enumeration, which counts only these bound- 
aries. As one can see, either upper index decreases or the lower one increases when going 
through the sequence. Now we introduce new pair of indices, which distinguish these two 
situation. Each group within which the lower index does not change, such that for some % we 
have N Pt _ x > N p . = N Pi+1 = ■■■ = N n+C ^ n y 1 = n > N Pi+c ^, is uniquely characterized by 
number n, 1 < n < N and cardinality c(n) G {0, . . . , m + 1}. This means that the particle 
number n appears c(n) times in the correlation function. It is convenient to introduce a 
pair of indices (n, a), where index n is the number of particles arriving at given boundary 
and index a, < a < n — 1, labels the position of given boundary within the group. Then, 
instead of the notation r/(A;) 6 Bj we use r^ n ' a ^ G B^ n ' a \ implying that for each transition 
between two A^-boundaries, in which the particle number does not change, the second in- 
dex a increases by 1, while in each transition within single iVp-boundary, which effectively 
reduces the number of fictitious particles by one, index n decreases by one. As a result, the 



r.h.s. of (III. 41) can be rewritten in a more uniform way 



const x det[*^_,(r^ u ')]i<fc,i<M 

Ni 

x n^^^^" 1 ' ^^^]!^^™ ( IIL55 ) 



71=1 

c(n) 

X 

a=l 



W det [J r (n,o),(n,a-l) (^l , A™ )] l<k,l<n- 



We are in position to apply Theorem 4.2 from [T7]. It states that the measure (III. 55) 
is determinantal and gives a recipe of construction of the correlation kernel for given ini- 
tial conditions. Specifically, let us define function 0( ra i' a i).(™2,a2) Q f transition between the 
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boundaries £ (ni ' ai) and B {n ^ a2) 

0(ni,oi),(n 3 ,a fl ) = JT^^^^ * ni+1 * J" (ni+ltC(ni+1))t(ni+1)0 ) (111.56) 

* ' ' ' * 4>n 2 * •^(n2,c{n2)),(n2,a 2 )-i 

where we used a definition of convolution 

(a*b)(x,y)= ^ a(x,z)b(z,y), (111.57) 



> (n,o) 
- 



with the summation in z performed over the points of the boundary B^ n,a \ which is between 
the boundaries where the indices x and y live, and 

y(n,a) = ^aUmfl)) + (nL5g) 

where ^tl? = ^nI-i( t{Ni ' 0) )- The argument r( n ' a ) of lives on #( n ' a ) due to the 

convolution with the function 0(( n ' a )>( Ari >°)). For the cases when c(n) = we formally define 

•7>,o),(n,o)0r,2/) = &x,y 

Consider matrix M with matrix elements 

M k>l = (fc * 0((M*)),W,o)) * (r ,- 1)5 (HI.59) 

where we can omit the dependence of t^ _1 (z) on the label i of the iV-boundary Bj. If the 



matrix M is invertible, the normalizing constant of the measure (III. 55 ) is equal to (det M) 1 . 



According to the Theorem 4.2 from [T7], the correlation kernel of (III. 55) is as follows 



K{b^ ni ' ai \ri)M n2,a%) {r 2 )) 

= E £ *fcf } (n) [M-% (4n * ^OW™,-))) ( r « T2 ) (IIL60) 
fc=l i=l 

_0(ni,ai),(n a ,o fl )( Ti)T2 ^ 

Furthermore, if the matrix M is upper triangular, the derivation of the kernel is significantly 
simplified. In this case we construct the set of functions {$[ n ' a ' ) }, which form a basis of the 
linear span of the set 

* 0(l lC (l))-(n,a) )(T O ; r ), . . . , (0 n * ^(n,c(n)),(n,a) )(r n- 1 ? r)}> (m . 61) 

fixed by orthogonality condition 

^r\r)^r\r) = <Sy, i,j = 0,...,n- 1. (111.62) 



rez 
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Then the kernel takes the following form 
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K(6^Vi)^ (n2 ' a2 V 2 )) = Y,^t 1 -f^t 2 -k\ r ^ - (p {ni ' aiUn2 ' a2) . (111.63) 



fc=i 

As a result we have: 



Proposition III. 9 Given densely packed initial conditions 

b° = (0,-1,..., -N +1) (111.64) 



the correlation kernel of the determinantal measure (III. 55) has the form 



21 l u l > °2 J J 27rii> J 2iriw (w-v)(l/v+l/-n 2 -l) 

Ti r ,„ (III. 65) 



-l(n 2 > ni) / 



cfco (l-p^i)«i-*aw«i— 2 



27ri«) 2 (W-I) n 2-™1(1/U) + 1/7T2-1)' 

To 

where b { ^ ai) = (x^U) E B [rH ^\ i = 1,2 and tt 2 = ti 13 * 2 (b k n \) . 

Proof. We first introduce function F n (b) defined by an integral representation, similar to 
the one of F n (b), with different integration contour. 

F n (x,t) = ^(f ^(q + lY^-wyv. (IH.66) 
27ri/ ro i w V wJ 

One can check that this function has the following properties: 

<)) n *F k = -F k+ i, (111.67) 

F{n,a),(n,a-l) * F k = F k , (III. 68) 

and 

F k (x,t) = 0, when k < and - ac > t (111.69) 

Note that the choice of the contour T^i ensures uniform convergence of convolution sums, 
which may extend to x = — oo and t = oo. Therefore one can interchange summation 
and integration, from where the formulas ( III.67[III.68 ) follow. The choice of the contour 



becomes relevant for F k with positive k as in this case there is a pole at w = 1, which must 
be placed inside the contour. One also must keep in mind that the convolution with <p n 
applied to the function of a point at B^' ^ results in a function of a point at B^ 1,0 \ while 
the convolution with J r ( n ,a),(n,a-i) yields the transition from B ( - n ' a ~ 1 ^ to B^'^. Since F k = F k 

36 



for k < 0, *S^-?( r ) = (-l) Wl_ ^-Jv 1+ z(& (Ari ' 0) (r) - 6?), and hence, using (|lIL67|lIL68|, we 



have 



(-ir~ l F_ n+l (b^ a \r)-b^). 



Then, the elements of the matrix M defined in (III. 59) are 

F_ k+l+l {b^\T k k ^)-b»). 



M, 



k,l 



(fe,c(fc))/ k-l\ 



k-l 



(111.70) 



(111.71) 



It follows from the definition of 1 and formula (III. 69 ) that M ky \ = when k > I and Mij 



1. Therefore the matrix M is invertible and upper triangular and we can straightforwardly 
go to the orthogonalization procedure. 



Substituting the initial conditions (III. 64) we obtain 

t (n,a) 



(n,a) 



r) 



1 

2rd 



dw ( 1 — p 



1 



(w - lfw 



k x (n,a) +n _ k _ 2 



(111.72) 



To, 



where t («.«)) = &(™. a )( r ). It is not a surprise that this is the same function, as the 

one obtained in the case of single iV-boundary. Its argument lives on single boundary B^ n ' a \ 
and the orthogonalization procedure referring to this boundary feels no difference with the 
previous subsection: 



(n,a) 



^ 1 La fi v ~ 1 
t = <t> dv \ 1 — p 

; 2vri / V v 



,(n,a) 



' (v - \)-i- l v i- n - x 
(l/7r s (6(r)) - l)v + 1 



(111.73) 



Apparently, the double integral part of the kernel coincides with the one obtained in previous 
subsection as well. We only need to derive an explicit expression 



end we note that we start the series of convolutions in (III. 56) with applying them either to 
^ r (n2,c(n 2 )),(n2,a 2 ) OT ) if c ( n 2) = 0, to 4> n , 2 . These functions can also be expressed in terms of 
F k (x,t). Specifically, the expression for <p n obtained in the previous subsection is 

Mn,r 2 ) = -7r(6(™ +1 ' c W)'°(r 2 ))A(& (n ' 0) (n) -& {n+1 ' c(n)) '°(r 2 )) 



(111.74) 



and from (III. 42) 



^(n 2 , C (n 2 )),(n 2 ,a 2 )(r 1 ,T 2 ) = n (b^>°)F (b^ (n) - b^>°(r 2 )). (111.75) 

Therefore we can use formulas ( III.67|III.68 ) for convolutions, which show that the lower 
index of the function F k increases by one and the function itself picks up a minus sign every 
time the number n decreases by one. Finally we have 



i ) (ni,ai),(n 2 ,a 2 ) _ (_l) n i~ n 2 ff ^( n2 ' a2 )'°)_p ^™2,c(n 2 ) _ ^(n 2 ,a 2 ),0^ 



(111.76) 
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which again coincides with the expression obtained in single iV-boundary case. As a result 



we arrive at the kernel expression (111.65). ■ 

Finally adopting the arguments from the end of the previous subsection for the collection 
of the boundaries 3j$ , . . . , B k ^ n m we arrive at the Fredholm determinant expression, stated 



in the theorem |II.4| For the sake of mathematical rigor one would have to analyze the 
convergence of the series obtained (i.e. the properties of the operator K). Similar analysis 
however has been done in many papers and we address the reader to them [3], [0 [8j EH]- 



IV. ASYMPTOTIC ANALYSIS OF THE CORRELATION KERNEL 



Now we use the parametrization of the space-time plane discussed in subsection II B 
Below we evaluate the scaling limit of the correlation kernel, suggesting that the arguments 
of the kernel are associated with a pair of boundaries and particle numbers fixed by choosing 



two points at the path (11.22 )-(II. 24) being at the distance of order of correlation length from 
each other. 



Lemma IV. 1 Let us fix two points at the path (11.22 )-(II. 24) in the 6 — v plane 



ri = r + UiL 1/3 , 



(IV. 1) 



where i — 1,2 and correspondingly set rij = [Lz/(r»)] and 6i = 6(ri). Let us consider two 



boundaries B 1 = B\ and B 2 = B\ which approximate smooth curves according to (11.18) with 



the parameters 81,62 fixed above. Then (11.19) define the curves approximated by boundaries 



B\ and B\ 2 corresponding to particles n\,n%, respectively. For the coordinates 7} of points 
on the boundary we also suggest the scaling 



r t /L = Xi = X(n) + SiL 2/3 , 



(IV.2) 



with Ui,Si fixed as L — )■ 00 and the function x( r ) defined in the subsection II B as a deter- 



ministic part of the random variable x, obtained as a solution of the equation (11.20) given 
6{r) andvir). Then 



lim L 1/3 K(b 1 ni (xi);b 2 n2 (x2)) ~ K f T 2 {-K 2 )K Aivy2 (K c u 1 , K f Si, k c u 2 , K f s 2 ) 



L— >oo 



(IV.3) 
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where in the r.h.s. we have the extended Airy kernel (11.29), k c and Kf are the model- 
dependent constants ( TLSUuLSi ) and 



T 2 (vr 2 ) 

The sign " 
minors. 



x 



u, vr B (6(r 2 )) = l 



(IV.4) 



(v^7(l + C (0 '%o)) + v^(l - C (0 '%o))) I tt b (6(t 2 )) = p 

means the equality up to the matrix conjugation, which does not affect matrix 



Proof. We introduce the following functions 

f{w;6,x) = ^ 6 ' X ] + X \n(q + p/w) 



+ u(r) In 



w — 1 



w 



C(e,x)-x 



h(w) = ln(l - l/w). 



(IV.5) 
(IV.6) 



To analyze the double integral part of the kernel Kq, we represent it as a sum 

oo 



fc=i 



where the functions ^Z™' , are given in ( III. 72 ) , ( III. 73 ) . Note that, instead of the index 



in the superscript characterizing the number of the boundary, we placed the notation for the 
boundary explicitly, to reflect the dependence of the functions on the form of this boundary 
and not of the others (here B means the first particle boundary, while the index n shows 
that we have to shift it n— 1 steps back in horizontal direction). In terms of above notations 
the integrals entering the summands become 

dw 



r ,i 27riw2 



-Lf(w;e 2 ,X2)-L 1 / 3 zh(w)) 



1 2vri (w-l)((l/ir 2 -l)v + l)' 



(IV.8) 
(IV.9) 



where z = kL 1//3 and 7r 2 = 7r B2 (6(r 2 )). To obtain the asymptotics of K , we first evaluate 



the integrals for ^"^(ti) and $^-fc( r 2) asymptotically as L — > oo and then perform the 
summation. 



Taking into account (IV. 2) one can approximate the function f(w; 9i, Xi) up to the terms 
of constant order by 

f(w; e^Xi) = f n M + L- 2/3 s t gM, (IV.10) 
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where we introduce the notations 



fr(w) = f(w; 6{r),x{r)) 



(IV. 11) 



and 



g r [w) 



d((r) 



+ 1 ) ln(g + p/w) + 



d((r) 



— 1 )lnw 



(IV.12) 



dx J \ dx 

where £(r) = ((9(r), x( r ))- The position of the double critical point of function f r (w), which 

satisfies ft(wo) = f"(w ) = is 



w (r) = 1 + 



2u(r) 



(IV.13) 



q(((r) -x(t))' 

Instead of the exponentiated functions we use their Taylor expansion at the points Wi = 
w (ri), with % = 1,2 for *"J_ fc (xi) and $J£_ fc (x 2 ) respectively. 



/r 4 H ~ / ri K) + ^fr (w )(w - W;) 

^r 4 H « 9n(wi) + g' ro (w )(w - Wi) 
h(w) ~ /i(iUj) + h'(wo)(w — tUj) 



(IV. 14) 
(IV.15) 
(IV.16) 



where in the coefficients of w-dependent terms we, without loss of accuracy, replace r\ and 
Wi by r and w = w (r ) respectively. We substitute these expansion into the integrals, and 
choose steep descent contours such that they approach the horizontal axis at the points w± 
and W2 at the angles ±7r/3 and ±27r/3 respectively. Changing the integration variables to 
£i = (w — w^L 1 / 3 f"'(wo)/2 we arrive at the integrals defining the Airy functions: 



Ai(a) 



°° eW3 dx 



t/3 



27Ti 



t exp 



X" 



— xa 



(IV. 17) 



As a result we have 



n\ — k 



n) 



exp [Lf ri (wi) + L 1/3 (si5( ri (wi) + zh(wi))) 



x Ai 



Wo 2 (L/;> )/2)V3 
zft/( w o) - s x cf (w y 



(IV.18) 



(/; n 'K)/2)V3 



( | ) n2,B 2 

ri2 — k 



(r 2 ) 



x 



exp (-Lf r2 (w 2 ) - L 1/3 (s 2 g r2 (w 2 ) + z/i(w 2 ))) 
K - 1)((1/tt»(6(x2)) - 1H + 1)(L/-K)/2)V3 



(IV.19) 
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The summation over k can be replaced by an integration over z. To perform the summations 
we use one more expansion: 



h(wi) = h{w ) - h'{w )w , {r )u i L- 1 ^ + 0(L- 2 / 3 ). 
Finally, taking into account that h'(wo) = \/{wq{wq — 1)), we obtain 

oo 



k=l 



X 



T 2 (7r 2 )L _1/3 fi; / e (L(/l/ i {u ' l)_/ ''2( M ' 2 )) +Ll/3 ( Sl 3( ,i 'i)- s 2s(«'2))) 

POO 

/ d\e XKc{u2 ~ Ul) Ai (A + KfSi) Ai (A + k } s 2 ) , 



where 



2V3 
2 1/ V(^o) 



and 



T 2 (tt 2 ) = - [g'(w Q )w ((l/n B (b(r 2 )) - l)w Q + I)]" 1 



Substituting 



"" " 2 7 g(v^7(l + C (0 ' 1) M) + v / ^(l-C (0 ' 1) (^o))) 

V (n>) (C (r ) - X (r„) C {0,1) (r )) (v^J - VtT 

-^'(r )C (1 ' 0) (r )(v^w-V7)] 
2g 3 7 5/2 



X 



£>o) 



#'(wo) 



g_V2 jym + C (0|1) M) + ygi - C ( °'%o))) 



(IV.20) 



(IV.21) 



(IV.22) 
(IV.23) 

(IV.24) 
(IV.25) 



(IV.26) 
(IV.27) 



1\fW> (%/P7 _ \/^) 
we have pI.30|II.31|IV4| ). 

Let us now evaluate the second part of the kernel given by the single integral, which can 
be written as 

dz exp [L(f ri (z) - fr 2 (z)) + L 1/3 (sig ri (z) - s 2 g r2 
2vri z((l/vr B (6(r)) - l)z + 1) 

The critical point of the exponentiated function is found to be 



z c = w = w(r ) 



(IV.28) 
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Then using the Taylor expansions we show that 



/riO) - fr 2 {w) « f n (Wl) ~ fr 2 {w 2 ) 



(IV.29) 



+ 
+ 
+ 



and 



6L 




r:' {^){u\-ui 


K(n>) 2 


2L 2 / 3 




O™o)(«2 ~ Ui 




2LV3 




si5f ri (wi) - S2g ra (w 2 ) 


+ (U 2 S 2 - UjSi) 


^- 1/3 4 



(2; - tt>o) 



2 - U7 J 



(IV.30) 



+ (si - s 2 )^ Q (w )(^ - w ) 

Substituting these expansions into the integral and integrating along the vertical line crossing 
the horizontal axis at wq we obtain: 



4k c (u 2 — «l) 



-K c Kf(s2U2 — Sl«l) 



(IV.31) 
(IV.32) 



y/ATTK c (u 2 - Ml) 

One can see that the first line of this expression exactly coincides with the factor before 



the integral in (IV. 21). Furthermore, its exponential part does not change the value of the 
determinants, so that it can be omitted. The second part can be rewritten using the formula 
from 130 



1 



'A-k{t' - r 

"OO 



3 -(5-€') 2 /4(r'-r)-(r'-r)(C+C')/2+(r'-r)3/12 



/oo 
-00 



(IV.33) 
(IV.34) 



where we should set r = k c Ui, t' = k c u 2) ^ = K/Si, £' = KfS 2 . As a result we obtain the Airy 
extended kernel 

L _1/3 T 2 (7r2)/t / i^Airy 2 (K C Ml, K/ S l5 K C U 2 , Kf8 % ) 



To finish the proof of the theorem II. 5 one has to prove the uniform convergence of the 



kernel in bounded sets and that the part of the sum (III. 34) coming from the complement 
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to these sets is negligible while the bound is uniform in L. For similar proofs we address 
the reader to J3J El |E1 [TT] . After that interchange of the sum and the limit is allowed. 
However we note that the limiting expression for the kernel still depends on which site 
b® 2 (72) is via the value of T2(vr2), which in turn depends on 7r s2 (6(r2)). To go from the sums 



(III. 34) to integrals we note that within every summation in an index r running through a 
boundary B the function Tg(7r B (6(r))) will enter linearly as a coefficient. It turns out that 
this coefficient amounts exactly to unit. This happens because the boundaries defined in 



(11.18) are locally straight and the r-dependent coefficient is averaged out on a smaller scale 
than the fluctuational one, which affects the resulting integral. The following lemma shows 
how the averaging works. After we apply it the statement of the theorem 2.2 follows. 

Lemma IV. 2 Suppose that b(r) e B, r = r + [sL 1/3 ] ; f }im (s) is a differentiable function 
and the following limit 

lim lM*f{b{r)) = T B (n)f lim (s) (IV.35) 

L— >oo 

holds uniformly in bounded sets s G [a, b] . Then, if the boundary is close to a continuous 
differentiable path in a sense (iLl^ILl^ , we have 



lim T = / fUs)ds (IV.36) 

r=ro+[aiV3] 

Proof. The proof is based on the fact that the order of the correction term accounting for the 
difference between the boundary on the lattice and its continuous differentiable counterpart 
allows one to consider the boundary as locally straight at the scales up to the fluctuation 
scale. This in particular means that in such a small scale, where the site-independent part of 
the limiting function can be considered as constant, the site-dependent part can be summed 
separately. It turns out that under this summation the site dependence exactly cancels with 
the slope dependence defined at the macroscopic scale, so that the remaining expression 
converges to integral of the site-independent part only. 

To be specific, let us divide the range of summation into bins of size eL 1 ^ 3 , where e is 
small, and perform the summation in two stages: first within each bin and second over all 
the bins. The first summation yields 

T +[e(n+l)L l /3] 

Mr)) * L- l '\NlT B {p) + NlT B (l))(fUen) + 0(e)), (IV.37) 

r^o + len.L 1 / 3 ] 
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where and N% and N% are the numbers of horizontal and vertical segments of the boundary 
within the summation range. Note that the fraction of these numbers, which corresponds 
to the slope of the boundary, being defined on the macroscopic scale persists up to the 
fluctuation scale, i.e. depends only on the value of x — h m L-s-oo tq/L: 

K = st^ etf/s ( 9«e,x)/dx + ^ (IV38) 

Nt = -Sx ^ -sit" r ^xV^X-n . (IV 39) 



From the explicit form of Yg(7r), (IV. 4), we have 

mmc±i Ta(p) _ s ^x)/9x-i Tem = x (IV40) 

i.e. (NyY(p) + A^|Y(1)) ~ eL 1 ^ 3 . Finally, after taking limit L — > oo, performing the second 
summation Z)i<n.<[(6-o)/e] £ fnm(. £n ) an d taking limit e — > we arrive at the desired result 

dTv36l). ■ 
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